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A Search for Vulcanoids with the STEREO Heliospheric Imager 

A. J. Steffi 1 , N. J. Cunningham 2 , A. B. Shinn 1 , D. D. Durda 1 , and S. A. Stern 3 

ABSTRACT 

Interior to the orbit of Mercury, between 0.07 and 0.21 AU, is a dynamically 
stable region where a population of asteroids, known as Vulcanoids, may reside. 
We present the results from our search for Vulcanoids using archival data from 
the Heliospheric Imager-1 (HI-1) instrument on NASA's two STEREO spacecraft. 
Four separate observers independently searched through images obtained from 
2008-12-10 to 2009-02-28. Roughly, all Vulcanoids with e < 0.15 and i < 15° 
will pass through the HI-1 field of view at least twice during this period. No 
Vulcanoids were detected. Based on the number of synthetic Vulcanoids added 
to the data that were detected, we derive a 3a upper limit (i.e. a confidence level 
> 0.997) that there are presently no Vulcanoids larger than 5.7 km in diameter, 
assuming an R-band albedo of pr=0.05 and a Mercury- like phase function. The 
present-day Vulcanoid population, if it exists at all, is likely a small remnant of 
the hypothetical primordial Vulcanoid population due to the combined effects of 
collisional evolution and subsequent radiative transport of collisional fragments. 
If we assume an extant Vulcanoid population with a collisional equilibrium dif- 
ferential size distribution with a power law index of -3.5, our limit implies that 
there are no more than 76 Vulcanoids larger than 1 km. 

Subject headings: Asteroids, Mercury, Image processing 

1. Introduction 

Interior to Mercury's orbit is a dynamically stable region where a population of small, 
asteroid-like bodies called Vulcanoids has long been hypothesized to exist, cf. review by 
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C ampins et al.l (ll996f ). This region, known as the Vulcanoid zone, extends from roughly 0.07 
AU to 0.21 AU (15-45 solar radii). Seen from Earth, objects in the Vulcanoid zone have 
maximum solar elongation angles of just 4-12°. The outer boundary of the Vulcanoid zone, 
at 0.21 AU, is set by dynamical instability. Objects with semi-major axes greater than this 
limit evolve onto Mercury-crossing or bits on 100-Myr timescales due to o rbital perturbations 
caused by both Mercury and Venus ( lEvans and Tabachnikl Il999l . 120021 ) . The inner edge of 
the Vulcanoid zone is less well-defined, but is set by the combination of the intense thermal 
environment and dynamical transport mechanisms such as Poynting-Robertson drag and 
the Yarkovsky effect. At 0.06 AU ( 13 solar radii) solar radiation is so in tense that even 



a pure iron body 10 km in diameter will evaporate in less than 4.5 Gyr f lLebofskyi 11975 



Campins et al.lll996l ). The time an object can survive against evaporation is a strong function 
of heliocentric distance, such that, at 0.07 AU, a pure iron body of just 2 km diameter 
will survive for the current age of the solar system. Poynting-Robertson drag extends the 
evaporation limit outward, as it can move a 2-km diameter object with p = 4 g cm -3 from 
0.08 to 007 AU in 4 Gyr, where it will evaporate (IStern and Durdal 120001 ) . 



The detection of one or more members of the putative Vulcanoid population is of interest 
as it would represent the discovery of a whole new class of solar system objects. If they existed 
at all, primordial Vulcanoids likely formed from the highest temperature condensates near the 
inner edge of the solar nebula, and they presumably would contain unique, highly refractory, 
chemical signatures as a result. In addition, a primordial Vulcanoid population might have 
affected Mercury's surface chronology. Based on the size-frequency distribution of craters 
on the Moon, Mars, Venus, and Mercury, it is thought that objects in the inner solar system 
were resurfaced du ring the period of the Late Heavy Bombardment (LHB) 3.9 Gyr ago 
( IStrom et al.ll2005l ). Vulcanoids removed fro m stable orbits in the Vulcanoid zone by non- 
gravitational forces like the Yarkovsky effect (jVokrouhlickyl Il999l ; IVokrouhlicky et al.l 120001 ) 
could have supplied a significant impactor population to M ercury after the LHB, making the 
surface appear older (ILeake et al.lll987l ; iHead et al.ll2007l ). 



1.1. Early Searches 



Le Verrierl ( 118591 ) first proposed that a small planet, or collection of planets, interior to 



Mercury could explain the observed precession of Mercury's orbit. This hypothesized planet 
was eventually given the name "Vulcan", and numerous searches were conducted in an 
attempt to find it. However, the proximity of intramercurial objects to the Sun makes them 
difficult to observe from the Earth. To overcome the observational challenge of looking for 
a faint object against a bright twilight sky, many of the early searches for the planet Vulcan 
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were conducted in the fleeting minutes of totality during solar eclipses. The first searches 
for Vulcan were done visually, with obvious li mitations . The first use of photographic plates 
to search for an intra-Mercurial planet was by iPerrind (119021 ). during the total solar eclipse 
of May 18, 1901, on the island of Sumatra. After analyzing the plates from this expedition, 
Perrine placed a limit of magnitude 5.0 (photographic) on the brightness of any planet 
interior to Mercury and concluded that, "...there are probably no bodies of appreciable size 
in the region close about the Sun, and that the cause of the disturbance in the motion of 
Mercury must be sought elsewhere." Subs equent observations during solar eclipses r esulted 
in an improved limiting magnitude of 8.0 ( jPerrine and Campbell! 119071 ; |Perrinelll909l ). 



In 1915, Einstein! (1191 5h showed that the precession of Mercury's orbit could be explained 
entirely by the then new theory of general relativity, thus eliminating the dynamical need for a 
massive planet Vulcan. However, the question of whether there are any small bodies interior 
to Me rcury's orbit remained unresolved. An archival search by ICampbell and Trumpler 
( 119231 ) using photographic plates obtained during the solar eclipse of 1922 for the purpose 
of measuring the deflection of starlight predicted by general relativity provided the tightest 
constraint of all the early searches: a photographic magnitude of 8.5. 



1.2. Modern Searches 
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More recently, ICourten et al.l (119761 ). summarizing a 10-year campaign to observe comets 
during solar eclipses, reported, "...data which indicate the possible existence of one or more 
relatively faint objects within twenty solar radii... and [ranging] from +9 to +7 in equiva- 
lent visual magnitude". Unfortunately, the nature of observing during total solar eclipses 
precluded any direct follow-up observations of these possible detections. Since these re- 
sults remain unpublished and have not been independently confirmed, it is difficult to assess 
whether the claimed de tections are real , and, if they are, whether the objects are Vulcanoids, 
sungrazing comets (cf. iBiesecker et al.l (120021 ) ). or some other type of body. 



In contrast to all other published Vulcanoid searches, iLeake et al.l (119871 ) conducted a 
search for Vulcanoids between 1979 and 1981 at 3.5 /im. By virtue of operating in the 
thermal infrared, this search was more sensitive to objects with low visual albedo. Leake et 
al. estimated a detection probability of 75% for an object with an L-band magnitude of 5, 
corresponding to an object diameter of 40-50 km. However, bad weather and the small field 
of view (FOV) of their instrument limited their search to a total of 5.8 deg 2 within 1° of the 
ecliptic-a small fraction of the Vulcanoid zone, as seen from Earth. 



69 



Prior to our work, the most complete search, in terms of depth and coverage, was 
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conducted by IDurda et al. (2000), using data from the LASCO C3 coronagraph on the 
SOHO spacecraft ( IBrueckner et al.lll995l ). which images a region from 0.02-0.14 AU (3.7- 
30 solar radii). Durda et al., examined a 40-day sequence of LASCO C3 images. Except 
for objects at the outer edge of the Vulcanoid zone with inclination >25°, all dynamically 
stable Vulcanoids should have passed through the instrument FOV during this period. No 
Vulcanoids were found, do wn to detection limit of V=8.0. For objects with a Mercury-like 
albedo and phase function (jVeverka et al.lll988l ). this limit corresponds to a diameter of 22 
and 65 km for objects at the inner and out er edges of the Vulcanoid zon e (or 36-106 km for 
an albedo of 0.05). Working independently, ISchumacher and Gayl (120011 ) also did not detect 
any Vulcanoids in LASCO C3 images down to a limiting magnitude of V=7. 



Subsequently, IDurda and Sternl (120031 ) conducted a Vulcanoid search using a visible 
wavelength imaging system flown aboard NASA F/A-18B a ircraft at an altitud e of 49,000 
feet. However, they were unable to improve upon the earlier IDurda et al.l (120001 ) results. 



Merline et al.l (120081 ) report preliminary results from their search for Vulcanoids using 



the Wide Angle Camera (WAC) of the MESSENGER spacecraft's Mercury Dual Imaging 
System (MDIS), while the spacecraft was in cruise to Mercury. Spacecraft pointing restric- 
tions limited observations to those with solar elongation >30°, i.e., the outer 45% of the 
Vulcanoid zone. No Vulcanoids were detected, down to a limiting magnitude of V=8, corre- 
sponding to an Vulcanoid diameter of 15 km. Subsequent analysis of these observations has 
brought this size limit down to 5 km, comparable to the results of our search (W. J. Merline, 
private communication, 2012). 



Finally, IZhao et al.l (120091 ) report a Vulcanoid search using 15 cm telescopes equipped 
with CCDs at two separate observatories in China during the 2008 total solar eclipse. They 
found "three unidentified star-like objects" in images from both telescopes, but the relative 
motion of these objects did not match that of a Vulcanoid. Both the angular size and 
sensitivity of this search are unclear, although they report that stars as faint as V=12.8 were 
detected. 



97 2. HI-1 Data and Processing 

98 For our search, we used archival data from the Heliospheric Imager (HI) instrument on 

99 NASAs STEREO spacecraft, available online from the STEREO Science Center ( |http : / / stereo-ssc . nasc 

100 The Solar TErrestrial RElations Observatory (STEREO) mission is designed to study coro- 

101 nal mass ejections (CMEs) from the Sun out to the orbit of the Earth and consists of two 

102 nearly identical spacecraft. STEREO-A orbits the Sun slightly interior to Earth's orbit, 



-5 - 



103 while the other sp acecraft, STEREO-B, orbits the Sun slightly exterior to Earth's orbit 

104 ( IKaiser et al.ll2008l ). Seen from the Sun, the angular separation between the Earth and the 

105 two spacecraft increases by 22.5° per year. Our search focused on data from STEREO-A, as 
we it has more stable pointing than STEREO-B. 

107 The HI instrument, itself part of the Sun Earth Connection Coronal and Heliospheric 

iob Investigation (SECCHI), consists of two separate imagers: HI-1, with a 20° square FOV 

109 centered in the ecliptic plane at a solar elongation o f 14°, and HI-2, w ith a 70° FOV centered 

no on the ecliptic plane at a solar elongation of 53° (lEyles et al.l 120091 ) . The HI instruments 

in use a series of baffles to prevent scattered light from the Sun from reaching the detector; as 

n2 a result, HI is significantly more sensitive to faint targets like CMEs (or Vulcanoids) than 

in traditional coronagraphic imagers. The HI-1 A imager is particularly well-suited to searching 

H4 for Vulcanoids, as all objects on dynamically stable orbits will be in the FOV at their eastern 

H5 elongation, provided the following inequality is satisfied: 



a\J 1 — e 2 cos % > r tan (3 (1) 

n6 where a, e, and % are the semi-major axis, eccentricity, and inclination of the Vulcanoid's 

ii7 orbit; r m 0.97 AU is the heliocentric distance of the STEREO-A spacecraft; and = 4° is 

us the angle between the sun and the inner edge of the FOV. This inequality is satisfied for all 

n9 orbits with a > 0.071 AU, e < 0.15, and i < 15°. 

120 The HI-1 detector is a 2048x2048 pixel CCD. To obtain sufficient signal-to-noise (S/N) 

121 for observing faint CMEs, HI-1 typically integrates for 1200s. A single 1200s exp osure would 



122 be sig nificantly degraded by cosmic ray hits, which affect roughly 45 pixels s 1 ( lEyles et al. 



123 120091 ). Instead, HI-1 typically acquires 30 separate 40-second exposures. Since the data vol- 

124 ume required to transmit all these images to Earth at full resolution would be prohibitive, the 

125 individual images are processed and combined onboard the spacecraft. The image processing 

126 includes a cosmic ray removal algorithm that compares each new image with the previous 

127 image on a pixel-by-pixel basis. If the value in a given pixel exceeds that in the previous 
i2B image by more than 5cr, where o is the predicted standard deviation based on the number 

129 of photoelectrons detected, that value is replaced with the value in the previous image. As 

130 seen from STEREO-A, the apparent motion of a Vulcanoid between 40s exposures is < 5". 

131 Given the 35"size of the HI-1 pixels, a full-width at half-m aximum value of the PSF of 3.34 



132 pixels in the X-direction and 2.96 pixels in the Y-direction (IBewsher et al.ll2010l ). the rate of 

133 apparent motion for Vulcanoids, and that the signal from the solar F-corona is much brighter 

134 than the signal expected from a Vulcanoid, this cosmic ray removal algorithm does not affect 

135 the ability of the instrument to detect potential Vulcanoids. 
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136 After the cosmic ray removal, the bias level, as determined from a column in the CCD 

137 underscan region, is subtracted from the data and the images are binned by a factor of two. 

138 The resulting sequence of 1024x1024 pixel images are then summed and compressed using 

139 a lossless algorithm, and in general, only these processed, summed images are sent back to 
wo Earth. In normal operations, these 1200s combined images are obtained with a cadence of 
i4i 40 minutes. 



142 2.1. Image Processing 

143 HI data are available in three levels of processing: L0, LI, and L2. L0 data consists 

144 of the raw, uncalibrated image data. LI data have been processed to correct for flat-field 

145 effects, image alignment, and the shutterless readout. As shown in panel A of Fig. [TJ both 
we L0 and LI data are dominated by light scattered from the solar F-corona (zodiacal light). 

147 In a typical 40x40 pixel region near the right-hand edge of the image, the median value is 

148 118 DN/s. The signal from the F-corona is fairly stable with time, so it can be effectively 

149 removed from the data by subtracting a composite image formed from data taken over a 

150 sufficiently long baseline. In the L2 data, the composite image to be subtracted is created by 

151 averaging over the lowest 25% of pixel values obtained during either a 1- or 11-day running 

152 window at each pixel in the image. The result of L2 processing can be seen in panel B of 

153 Fig. [TJ The same 40x40 pixel region has a median value of 0.2 DN/s in the L2 data. We 

154 used L2 data from the HI-1 instrument of STEREO-A (henceforth HI-1A) with the 11-day 

155 running window in our search. 

156 With the F-corona effectively removed, CMEs and coronal streamers are readily appar- 

157 ent in the HI-1 data. Although the STEREO mission is designed to study solar phenomena 

158 such as these, for the purpose of searching for Vulcanoids they constitute an additional source 

159 of background "noise" . If not removed, this additional background, which is more significant 
wo at smaller angular distances from the Sun, will result in reduced sensitivity to Vulcanoids 
lei with small semi-major axes. We attempted to remove this two-dimensional, time-varying 

162 background from the L2 data by subtracting off the "sky" level in a circular annulus centered 

163 on each pixel. Based on an average PSF FWHM value of 1.57 pixels (binned), we adopted 

164 radii of 4.7 and 7.9 pixels (3x and 5x the FWHM) for our sky background annulus. This 

165 method removes signal from diffuse features while preserving point sources. A sky-subtracted 
lee image is shown in panel C of Fig. [TJ After this sky-subtraction, the median value in the 

167 40x40 region fell to 8xl0~ 3 DN/s. 

168 Since the pointing of the STEREO spacecraft was held fixed relative to the Sun, back- 

169 ground stars appear to drift through the HI-1 FOV. To facilitate the identification of objects 
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170 moving relative to the inertial frame, we co-registered all sky-subtracted L2 data obtained 

171 during a two-day window (typically a series of 72 images). The first image in the sequence 

172 was used to define an astrometric reference frame. We then calculated centroid positions 

173 for several hundred stars common to both images and used these to derive a third degree, 

174 two-dimensional polynomial transform that maps the second image into the frame of the 

175 reference image. This polynomial technique handles optical distortions near the edges of 

176 the instrument FOV better than a simple linear shift. After the data were co-registered, we 

177 subtracted a 9-image running median from each image to reduce the signal from any sources 

178 that remained fixed in inertial coordinates, e.g., stars. An example of the final processed 

179 image is shown in panel D of Fig [TJ 



wo 2.2. Synthetic Vulcanoids 

181 To quantify the detection efficiency of our search, we generated a population of 101 

182 synthetic Vulcanoids and added them to the L2 HI-1A data, with the appropriate level of 

183 Poisson noise, before our additional image processing. The orbital properties of the synthetic 

184 Vulcanoids were chosen to mimic those of the putative Vulcanoid population. During the 

185 search, the observers knew that synthetic Vulcanoids were present in the data, but they 
we did not know where the synthetic objects were, how many objects were in a given image 
187 (between 23 and 42 with an average of 33), or whether a given object that was detected was 
las real or synthetic. 

189 Each of the 101 synthetic Vulcanoids was randomly assigned a value for the longitude of 

190 its ascending node, argument of periapsis, and mean anomaly at epoch. Given the FOV of 

191 the HI-1A instrument, we limited the range of orbital eccentriciti es to 0-0.15 and inclinat ions 



192 to 0-15°, with the value for each chosen randomly. Prior work by lStern and Durdal (120001 ) has 

193 also shown that Vulcanoids are more likely to be found on highly circular orbits near the outer 

194 edge of the Vulcanoid zone. For the synthetic population as a whole, (e)=0.07 and (i)—7.7°. 

195 Half of the synthetic Vulcanoids were given a semi-major axis randomly-distributed between 
we 0.07-0.14 A.U.; the other half were given semi-major axes distributed randomly between 

197 0.14-0.21 A.U. For each image to be examined, we calculated the positions of the synthetic 

198 objects in detector coordinates, and if the object was within the instrument field of view, we 

199 added a 2-D Gaussian PSF (FWHM of 1.67 pixels in the X-dimension and 1.48 pixels in the 

200 Y-dimension) with Poisson noise to the image. The positions of the synthetic Vulcanoids, in 

201 detector coordinates, over the 40-day search period are shown in Fig [2j 



202 The apparent brightness of some Vulcanoids can change by nearly three orders of mag- 

203 nitude as they move through the HI-1A field of view. This is primarily a result of the large 
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204 range of phase angles at which they can be observed (roughly 20°- 160° for a Vulcanoids in 

205 the outer Vulcanoid zone and 60°- 120° for Vulcanoids in the inner Vulcanoid zone) as well 

206 as potentially significant changes in the object's heliocentric distance and distance to the 

207 STEREO-A spacecraft while they are in the field of view. The synthetic Vulcanoids were 

208 assumed to be spherical and within each group (a < 0.14 AU and a > 0.14 AU), they were 

209 assigned diameters that uniformly spanned the range of 0.7-50 km, in logarithmically-spaced 

210 steps. The total signal (in DN/s) from each synthetic Vulcanoid was calculated according 

211 to the following equation: 



[iAd 2 
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16r^r 2 A 2 G 



/ (/^) S ®( X ) R ( X > a ) H ( X ) dX 



(2) 



where A is the collection area of HI-1A (A = Tcdf/4 where di = 1.59 cm), d is the diameter of 
the Vulcanoid in km, ro = 1.496 x 10 8 km is the astronomical unit (AU), r is the heliocentric 
distance of the Vulcanoid in AU, A is the distance from STEREO-A to the Vulcanoid in AU, 
G is the system gain (15 photoelectrons / DN), X n /hc is the number of photons of wavelength 
A n per erg, £0 is the solar flux at 1 AU in erg s _1 cm -2 nm" 1 , R is the absolute reflectance 
of the Vulcanoid at X n and phase angle a, and H(X n ) is the instrument response (which 
includes the CCD quantum efficiency). Finally, fi = 0.93 is the correction factor for HI-1A 



betw een the theoretically predicted counts and counts actually measured ( IBewsher et al. 



20101 ) . Eq. [2] is integrated between 300-1100 nm, the range over which HI-1 is sensitive. The 
instrument response curve of HI-1 A, H(X n ), contains a broad central peak from 600-750 nm, 
roughly similar to the Johnson R filter. However, there are significant side bands at 400 nm 
and 950 nm (cf. Fig. 6 in Bewsher et al., 2010). For an object with a Mercury-like spectrum, 
roughly 90% of the detected flux will be in the central pass band, 4% in the 400 nm side 
band, and 6% in the 950 nm side band. 



The solar spectrum, S Q , was obtained from the HST CALSPEC database (jBohlin 



20071 ). The synthetic Vulcanoids were assumed to have spectral reflectance properties iden- 
tical t o that of the planet Mer cury, except with 1/3 the albedo, i.e., scaled to p# = 
0.047 (IWarell and Bergfora 120081 ) . This albedo is comparable to that of C-type asteroids 
( Tholen and Barucci 19891 ) . For objects with a different albedo, d oc p^ 2 . The absolute 
reflectance of Mercury, as measured by the MDIS and MASCS instruments on the MESSEN- 
GER spacecraf t at 11 wavelengths between 428-1013 nm at 5 phase angles from 35°-129° by 
Holsclaw et al.l (120101 ) is shown in Fig. [3J For wavelengths between 300-428 nm and 1013- 
1100 nm, we linearly extrapolate the MDIS/MASCS data. For reflectances at phase angles 
between 35°-129°, we interpolate the MDIS/MASCS data, while for phase angles between 
20 o -35 o /129 o -16 0° we scale the reflect ance spectrum at 35°/129° by the V-band Mercury 
phase function of iMallama et al.l (120021 ). Although we have used Mercury's reflectance spec- 
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238 trum in our analysis, our results are insensitive to the choice of spectrum for reasonable solar 

239 system alternatives, such as typical S- or C-type asteroids, when that spectrum is scaled to 

240 an R-band albedo of 0.047. 



241 3. Search Technique 

242 After processing the data and adding synthetic Vulcanoids as described above, we cre- 

243 ated movies from the images obtained during a two-day period. We then looped this movie 

244 back and forth, at variable frame rates, while visually searching for moving objects. When 

245 an object was found, the observer marked its position in one or more frames of the movie. 

246 This technique takes advantage of the ability of the human eye to detect motion; very often 

247 an object that was easily identifiable while cycling through the movie would be virtually in- 

248 distinguishable from background noise when looking at an individual image. After searching 

249 through all of the movies, the list of marked objects for each movie was compared to the 

250 positions of the synthetic Vulcanoids and main belt asteroids in that movie. Even with a 

251 population of just 101 synthetic Vulcanoids, with 70" pixels, there were multiple instances 

252 in which there were two synthetic Vulcanoids within a five-pixel radius of a marked posi- 

253 tion. These were automatically flagged for further review. Unless it was clear that both 

254 objects were marked in the same frame, only the brighter of the two synthetic Vulcanoids 

255 was recorded as being detected. Any marked position not matching the location of a syn- 

256 thetic Vulcanoid or main belt asteroid was flagged for further inspection. With one exception 

257 (discussed in Section I4T3"]) . these were subsequently determined to be false positives-usually 

258 the result of a interpolating a single "warm" pixel or a spurious grouping of random noise. 

259 Although there are now several years of STEREO HI-1 data, the labor-intensive nature 

260 of our search technique placed practical limitations on the amount of data that could be ana- 

261 lyzed. The longest synodic period between an object in a dynamically stable Vulcanoid orbit 

262 (i.e. having a semi-major axis < 0.21 AU) and an observer at 1 AU is 39 days. We therefore 

263 selected a 40-day period from 2008-12-10 to 2009-01-19 to examine in our search, based on 

264 the low solar activity and avoidance of the galactic center. (We had previ ously examined th e 



265 40-day period from 2009-01-19 to 2009-03-01 and found no Vulcanoids flSteffl et al.ll2010h . 

266 but we used a less effective sky subtraction technique resulting in a less sensitive search than 

267 the one described in this paper). For any given 40-day period, a Vulcanoid at the outer edge 

268 of the Vulcanoid zone will spend roughly 15 days continuously in the HI-1 A field of view 

269 during its single pass through elongation. In contrast, a Vulcanoid at the inner edge of the 

270 Vulcanoid zone will spend roughly 30 hours in the HI-1 A FOV per elongation, but will pass 

271 through elongation 5 or 6 times. 
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272 As has been pointed out by many others before us, given a careful observer, the human 

273 eye is often better at finding low S/N point sources with fewer false positives than automated 

274 search algorithms. This is especially true given the relatively dense stellar backgrounds in the 

275 HI-1 images and the large range of apparent motions that Vulcanoids can exhibit (anywhere 

276 from no apparent motion at all up to a few hundred arcseconds per hour, positive or negative, 

277 in both right ascension and declination). However, using human observers to manually look 

278 for Vulcanoids introduces a certain stochastic element to the search. To minimize the chance 

279 that a faint Vulcanoid might be missed due to observer fatigue or random chance, each movie 

280 was independently examined by four separate observers (AJS, NJC, ABS, and DDD). 

281 Finally, we note that based on this wo rk and previous searches, e.g. the search for 



282 small satellites of Pluto of ISteffl et al.l ( 120061 ). adding synthetic objects to the data not only 



283 provides a quantitative means of estimating the detection efficiency but actually results in 

284 a more sensitive search. The tedium of searching through large numbers of images for point 

285 sources marginally brighter than the background while failing to find any can quickly lead 

286 to observer fatigue and a less thorough search. However, the simple positive reinforcement 

287 of finding an object in the data, whether synthetic or not, can help keep a human observer 

288 focused on the search, despite its repetitive nature. 



4. Results 
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No Vulcanoids were detected in our search. Although we did not detect any Vulcanoids, 
we observed a variety of other solar system objects: the planets Mercury, Venus, Uranus, and 
Neptune ; the comet 67P/Chu ryumov-Gerasimenko; more than 30 Kruetz-family sungrazing 
comets (IBiesecker et al.l 120021 ) ; the non-group sungrazing comets C/2008 Y12 (SOHO) and 
C/2009 Al (STEREO); and numerous main-bel t asteroids, some as faint as V=13.8, as 
determined by the JPL HORIZONS ephemeris (jGiorgini et al.lll996l ). In addition, 58 of 
the 101 synthetic Vulcanoids added to the data were detected by at least one of the four 
observers, which we use to place upper limits on the size of any putative Vulcanoids. Finally, 
as discussed below, we initially classified the sun-grazing comet C/2008 Y12 (SOHO) as a 
candidate Vulcanoid. This resulted in an additional search through all images from the 
HI-1 instruments on both STEREO-A and STEREO-B over the period from 2008-12-08 to 
2009-03-01. We found no Vulcanoids in this additional search but cannot place quantitative 
limits on this additional non-detection, as no synthetic Vulcanoids were added. 
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4.1. Insensitivity to Orbital Parameters 

The detection efficiency of Vulcanoids, whether real or synthetic, will clearly be a strong 
function of the object's size, as the larger the object, the brighter it appears. It is also 
possible that our search might be biased towards detecting objects with a certain range of 
orbital parameter values. For example, in our prior Vulcanoid search we found that given 
two synthetic Vulcanoids of the same size, the o ne with the greater orbital semi-major axis 



309 was more likely to be detected (jSteffl et al.ll2010l ). To identify whether any such biases were 



310 present in our current search, we plot the orbital parameters of the synthetic Vulcanoid 

3n population as a function of diameter in Figs. H] and Synthetic Vulcanoids that were 

312 detected by at least one observer have filled symbols and those that were not detected at all 

313 have open symbols. From these figures, it is immediately clear that there is an extremely 

314 strong dependence on object size, with all 57 synthetic Vulcanoids larger than 5 km in 

315 diameter being detected by at least one observer and only 1 of the 44 synthetic Vulcanoids 

316 smaller than 5 km in diameter being detected (and then, only by one observer). However, 

317 there is no obvious correlation between whether an object was detected and any of its orbital 
3ib parameters. 

319 To quantitatively verify our search's insensitivity to a Vulcanoid's orbital parameters, we 

320 divided the synthetic Vulcanoids into two populations based on whether they were detected 

321 or not and conducted a series of statistical tests. These are described in more detail in the 

322 Appendix. None of these tests showed any significant deviations from a random uniform dis- 

323 tribution or correlations between pairs of parameters. This is a somewhat surprising result. 



324 Evidently the sky annulus subtraction technique described in Section I2.1[ has removed the 

325 bias towards detecting objects with larger semi-major axes observed in our initial search. 

326 Given the relatively low number of synthetic Vulcanoids (101), it was not possible to in- 

327 vestigate whether any statistically significant groupings of three or more orbital parameters 
32B exist. 



329 4.2. Detection Efficiency vs. Object Size 

330 We present the number of synthetic Vulcanoid detections for each observer, X obs , as a 

331 function of size in Table [TJ We have grouped the objects into 13 size bins, corresponding to a 

332 factor of \/z in diameter. As might be expected, there is some disparity between the detection 

333 rates of the four observers, with observer 1 detecting 58 of the 101 synthetic Vulcanoids 

334 and observer 4 detecting 50 of the 101 synthetic Vulcanoids. Surprisingly, though, all the 

335 differences between the observers occurred for objects 4-8 km diameter. All four observers 

336 detected all 45 objects larger than 8 km and none of the 38 objects smaller than 4 km. For 
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337 each observer, the detection efficiency, e, in a given size bin is simply X f, s /N. We plot e for 

33B each of the four individual observers as well as e avg , the average detection efficiency of all 

339 observers vs. object diameter in Fig. EJ lcr confidence intervals for the e avg are derived from 

340 Monte Carlo analysis, as described below. We find e avg = O.SOOIq^ (32/40) for objects 

341 between 5.7-8.0 km diameter and e avg = 0.125lo!o5i (4/32) for objects between 4.0-5.7 km. 

342 Given the individual detection efficiencies, we calculate the probability, Pdet, that at least 

343 one of the four observers would detect a single object, as a function of size. Subject to the 

344 assumptions stated in Section I2T21 (namely a R-band Vulcanoid albedo of p#=0.05) we find 

345 that our search should have detected all Vulcanoids larger than 5.7 km in diameter. However 

346 this approach ignores the statistical uncertainty in the observed values of e, which arises from 

347 the small number of synthetic Vulcanoids in each size bin. 

348 We use a Monte Carlo model to quantitatively determine the effect of statistical uncer- 

349 tainty on our limits. For each size bin and observer (or the sum of the four observers, in the 

350 case of e avg ) we set the detection efficiency, e, equal to a pseudo-random number generated 

351 over the interval [0, 1], inclusive. Using the probability mass function for the binomial dis- 

352 tribution, we calculated the probability of detecting exactly X obs synthetic Vulcanoids, given 

353 N Vulcanoids in the size bin and e. We then generated a second pseudo-random number 

354 between and 1, and if it was less than or equal to this probability the test was considered 

355 successful at reproducing the observations and we recorded the value of e. This process 

356 was repeated until 10 6 values for e were recorded, la confidence intervals were produced by 

357 sorting the recorded values of e into ascending order and finding the minimum contiguous 
35B interval (ehigh ~ e iow) that contained 68.3% of the e values. These confidence intervals are 

359 shown in Table [TJ 

360 To determine the probability that a single Vulcanoid would be detected by at least one 

361 observer, we randomly selected one of the recorded values for e for each observer and tested 

362 whether this value was less than or equal to a pseudo-random number, i.e. if that observer 

363 failed to detect the object. If this test was true for all four observers, the Vulcanoid escaped 

364 detection. We repeated this 10 6 times and counted the total number of times the Vulcanoid 

365 was detected. The probability of a single Vulcanoid being detected by least one observer, 

366 given statistical uncertainties, P r det , is shown in Table [TJ We note that the values for P' det 

367 are significantly higher than Pdet fc> r Vulcanoids below 5.7 km in diameter. This is due to 
36B the fact that, given the small values for N, even if X G fc s = 0, it is not possible to exclude 

369 a non-zero value for e with a high degree of confidence. For example, observer 4 detected 

370 of 8 objects between 4.0-5.7 km, but even a value of e = 0.31 would yield this result 5% 

371 of the time. Thus, when at least one observer has an observed detection efficiency, e, close 

372 to zero, P' det is likely too high. The converse is also true: when at least one observer has 

373 e close to unity, P' det is likely an underestimate. With this in mind, we are able to place 
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374 a 3cr upper limit (i.e., a confidence level >0.997) on the existence of any Vulcanoids larger 

375 than 5.7 km in diameter. As stated above, our size limits are dependent on the value of the 

376 R-band albedo, p# = 0.05, we assumed for Vulcanoids. If the actual albedo of Vulcanoids is 

— 1/2 

377 different from this value, our diameter limits scale as p R . 



4.3. C/2008 Y12 (SOHO) 

379 In our search, we found one object that did not correspond to a synthetic Vulcanoid 

380 or any known object in the IAU Minor Planet Centers database (Gareth Williams, private 

381 communication, 2011). This object was first detected in an image from 2008-12-21T17:44 UT 

382 and appeared in 15 consecutive HI-1A images before its apparent motion carried it outside 

383 the instrument's FOV. Although the object's motion was similar to several to the synthetic 

384 Vulcanoids, it appeared somewhat diffuse and brightened as it moved closer to the sun. We 

385 thus suspected that this candidate Vulcanoid might be a previously undetected sungrazing 

386 comet. However, given the short 9.3-hour arc and 70" pixels it was not possible to obtain 

387 a unique orbital solution. We therefore began a concerted effort to see if the object could 

388 be re-acquired in images from the ffl-1 instruments on either STEREO-A or STEREO-B 

389 obtained over the 83-day period from 2008-12-08 to 2009-03-01. These data were processed 

390 and searched using the techniques described above, but without the addition of any synthetic 

391 Vulcanoids. Ultimately, we were not able to reacquire the object and no other candidate 

392 Vulcanoids were detected. Since we were specifically looking to reacquire the candidate 

393 object in the data, this secondary search was likely not as sensitive to Vulcanoids in general 

394 as our primary search. However, lacking any synthetic Vulcanoids to detect, we can not 

395 quantitatively verify this hypothesis. Given that in our primary search, all four observers 

396 detected all 45 synthetic objects larger than 8 km, it seems likely that were such an object 

397 present in this extended data set, it would have been detected. This gives us additional 

398 confidence that some unknown systematic effect unique to either the HI-1A or the time 

399 period from 2008-12-10 to 2009-01-19 is not preventing the detection of a real Vulcanoid. 

400 Subsequently, we learned of the comet C/2008 Y12 (SOHO), which was discovered in 

401 7 images from the LASCO C2 coronagraph on the SOHO spacecraft, starting on 2008-12- 

402 22T16:54-less than 14 hours after the Vulcanoid candidate was last seen in the HI-1A images. 

403 Th e predicted posi tions of C/2008 Y12 (SOHO) in the HI-1A images using the initial orbital 



fit ( iMarsdenl 120091 ) were not a good match to our candidate object. However, we were able 
to find a single set of orbital parameters that yielded average positional residuals of 37" for 
the Vulcanoid candidate in the STEREO HI-1A images and 3" for C/2008 Y12 (SOHO) in 
the SOHO LASCO-C2 images-less than half the size of the respective detector pixels. After 
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408 reporting these results to the IAU Minor Planet Center, a new orbital fit was published 



Williams! ( 120121 ). However, this published fit assumed an eccentricity of 1.0, and as a result, 



it is significantly less accurate than our orbital fit shown in Table [2J 



Discussion 
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Although it is dynamically stable, the Vulcanoid zone is in a rough neighborhood. 
It occupies a comparative ly small volume o f sp ace and due to it s prox imity to the sun, 
orbital velocities are large. iLeake et all (119871 ) and lStern and Durdal (120001 ) showed that even 
for a wide variety of assumptions about the size and material properties of a hypothetical 
population, Vulcanoids will experience significant collisional disruption and erosion unless 
the mean orbital ecc entricity of the population is extremely small (~ 10~ 3 ). However, 
dynamical studies by Evans and Tabachnikl (120021 ) showed that test particles between 0.09 
and 0.20 AU on initially circular orbits with zero inclination evolved dynamically such that, 
after just 100 Myr, the population had a mean eccentricity of (e) = 0.0935 and a mean 
maximum eccentricity over the 100 Myr of (e max ) = 0.15. This increase in mean eccentricity 
is due to the interaction of various orbital resonances with Mercury as its orbit evolved 
over the simulation. For a Vulcanoid population with mean eccentricity of (e) = 0.0935, 
typical encounter velocities between Vulcanoids will be in the range of 10-20 km s _1 and 
collisional disruption/erosion is extremely efficient. For an initial Vulcanoid population with 
a power-law index of -2.5, 300 objects larger than 1 km radius, and (e) = 0.1024, high- 
energy colli sions with small bodies reduced all objects to debris, no larger than 1 km radius, 
in 1.2 Gyr (IStern and Durdal l2000l ). For an initial population of 10 4 objects with r> 1 km, 
this collisional destruction timescale was just 1.6 Myr. 

However, radiative transport mechanisms like the Yarko ysky effect would have qui ckly 
removed the smaller collisional remnants. In simulations by IVokrouhlicky et al.l (120001 ) all 
Vulcanoids 1 km in diameter and smaller were removed from the Vulcanoid zone by Yarkovsky 
drift, with typical lifetimes of 100 Myr for 0.1 km objects and 500 Myr for 1 km objects. 
Even objects up to 10 km in size showed strong depletion (80%-100%) over the age of the 
solar system. By removing the small body population, the Yarkovsky effect might greatly 
reduce the efficiency of collisional erosion in the Vulcanoid zone. To date, there have been 
no studies combining the effects of collisional evolution and Yarkovsky drift in the Vulcanoid 



zone. 



439 If we assume that the present-day Vulcanoid population bot h exists and obeys a steady- 

440 state collisional fragmentation size distribution, dn = C*r~ S 5 dr JPohnanvl 19691 : Williams and Wetherilll 



19941 ). then our 3a upper limit implies that no more than 76 primordial Vulcanoids remain in 
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the Vulcanoid zone-the last remnants of a much larger initial population. Such Vulcanoids 
would all be larger than 1 km in diameter, as primordial objects smaller than this will have 
been removed from the Vulcanoid zone by Yarkovsky drift. However, since the combination 
of collisional evolution and Yarkovsky drift could result in either the complete depletion of 
a primordial Vulcanoid population or something less than that, our upper limits cannot be 
extrapolated to place any meaningful constraints on a hypothetical primordial Vulcanoid 
population. 

Finally, we note that if a Vulcanoid is subsequently discovered (hopefully below our de- 
tection limits) it is possible that it may be a scattered Near-Earth Object (NEO) rather 
than a primordia l Vulc anoid. Recent modeling of the orbital distribution of NEOs by 



Greenstreet et al.l (120121 ) suggests that up to 0.006% of the NEO population might end 



up in orbits with an aphelion interior to Mercury's aphelion. Given that limit, an NEO 
interloper might have an absolute magnitude as low as H ~ 21, corresponding to a diameter 
between 170-370 m, more than an order of magnitude smaller than the detection limits of 
our search. 
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460 be C/2008 Y12 (SOHO), and K. Battams for discussions about C/2008 Y12 (SOHO) and 

461 its orbit. Support for this work was provided by NASA Planetary Geology and Geophysics 

462 program through grant NNX09AD65G. 



463 A. Statistical Tests 

464 To determine if there were any statistically significant deviations in the detected and/or 

465 non-detected populations of synthetic Vulcanoids, we conducted a series of statistical tests. 

466 The first test was to divide the range of possible values for each orbital parameter into four, 

467 equally-spaced bins. If the parameter values are randomly distributed, then the expectation 

468 is that 1/4 of the total population should have parameter values in the range covered by 

469 each bin. The cumulative binomial distribution can be used to determine the significance 

470 of any deviation from this expected value. The largest deviation from a random uniform 

471 distribution occurred in the non-detected population: the fourth semi-major axis bin had 17 

472 objects, compared to an expectation value of 10.75-an event with only a 2.5% probability. 

473 While this may seem marginally significant, given that we conducted 24 separate tests, we 

474 would expect one of them to have as large a deviation nearly 46% of the time. Thus, we 
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475 conclude there were no statistically significant deviations from a random normal distribution. 

476 We also performed one- and two-sample Kolmogorov-Smirnoff (K-S) tests on the dis- 

477 tributions of orbital parameters for both the detected and undetected populations. A one- 

478 sample K-S test compares the empirical distribution function of a sample population with 

479 the cumulative distribution function of a reference distribution (in this case, a random uni- 

480 form distribution) and evaluates the significance of any deviations. None of the empirical 

481 distribution functions of the six orbital parameters for either of the two populations was 

482 significantly different from a random uniform distribution. Similarly, a two-sample K-S test 

483 compares the empirical distribution functions of two sample populations to determine if they 

484 are consistent with being derived from the same parent distribution. We found no statisti- 

485 cally significant difference between the parent distributions of the detected and non-detected 

486 populations for any of the six orbital parameters. 

487 Finally, we examined each of the 21 combination of the six orbital parameters and 

488 Vulcanoid diameter for both possible correlations. Both the Pearson linear correlation coef- 

489 ficient (which measures how well two variables are linearly related) and the Spearman rank 

490 correlation coefficient (which measures how well the relationship between the two variables 

491 can be described using a monotonic function) exhibited near-zero values for all parameter 

492 combinations. 
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Table 1. Synthetic Objects Detected by Observer 
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a Diameter assumes a Vulcanoid R-band albedo of 0.05. For other assumed values 

— 1/2 

of the albedo, doc p^, 

b Detection efficiency, averaged over the four observers, with la (68.3%) confi- 
dence intervals. 

Probability of at least one of the four observers detecting a single object, given 
statistical uncertainties associated with X fe s detections of N objects for each ob- 
server. 

d When X o ft s =0 for at least one observer, P^ e ^ is likely to significantly overstate 
the detection probability. This results from the fact that given the small value of 
N, a non-zero value for e j, s cannot be excluded with a high degree of confidence, 
even though in all likelihood, e ~ 



Table 2. Orbital Elements for C/2008 Y12 (SOHO) 



Orbit Solution 


q(AU) 


e 


i(°) 




con 


T (UTC) 


Marsden (2009) 


0.0533 


1.0 


154.34 


35.03 


140.06 


2008 Dec. 22.65 


Williams C2012) 


0.0658983 


1.0 


23.39475 


312.10406 


145.62274 


2008 Dec. 22.60543 


This work 


0.0659134 


0.970428 


23.38917 


312.10392 


146.56106 


2008 Dec. 22.62134 
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Fig. 1. — STEREO HI-1A images at various levels of processing using the same stretch. Panel 
A shows an image after Level 1 processing. The solar F-corona (zodiacal light) dominates 
the image. Panel B shows the same image after Level 2 processing to remove the solar F- 
corona. Streamers and other features of solar origin are clearly visible. All point-like sources 
are stars. Panel C shows the Level 2 image in panel B after our additional sky subtraction; 
the solar features have been almost completely removed. Panel D shows the image in panel 
C after subtracting the median of the aligned images to suppress sources that are fixed in 
inertial coordinates (e.g. stars). See the text in Section [27X1 for more details. 
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Apparent Positions of 101 Synthetic Vulcanoids 




200 400 600 800 1000 
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Fig. 2. — Positions of the synthetic Vulcanoid population in HI-1A detector coordinates 
during the 40-day search period. Shades of gray indicate the instantaneous apparent R 
magnitude of the synthetic object. 
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Reflectance data from Holsclaw et al., 2010 
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Fig. 3. — Absolute disk-integrated reflectance of Mercury at fi ve phase angles, as m easured 
by the MASCS and MDIS instruments on MESSENEGER Jrlolsclaw et all boioh . The 
MASCS/MDIS data have been linearly extrapolated to cover the full wavelength range over 
which the STEREO HI-1 instrument is sensitive. For synthetic Vulcanoids at phase angles 
outside the range (35 -129°) cov e red b y the observations, we scale the data by the V-band 
phase function of iMallama et al.l ( 120021 ). 
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Fig. 4. — Orbital inclination, eccentricity and semi-major axis of synthetic Vulcanoids, as a 
function of the object's diameter. Filled symbols represent Vulcanoids that were detected 
by at least one observer in our search, while open symbols were not detected by any of 
the four observers. There is no obvious correlation between a synthetic Vulcanoid's orbital 
parameters (a, e, and i) and whether it was detected. In contrast, all objects larger than 
5 km were detected, while only 1 object smaller than 5 km was detected. 
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Fig. 5. — Mean anomaly, Longitude of the ascending node, and argument of periheliob 
of synthetic Vulcanoids, as a function of the object's diameter. Filled symbols represent 
Vulcanoids that were detected by at least one observer in our search, while open symbols 
were not detected by any of the four observers. There is no obvious correlation between a 
synthetic Vulcanoid's orbital parameters (u, Q, and M) and whether it was detected. In 
contrast, all objects larger than 5 km were detected, while only 1 object smaller than 5 km 
was detected. 



-27- 



1.0 

^ 0.8 

o 
c 

CD 

'§ 0.6 

LU 

I 0.4 

o 

Q) 

Q 





X 


e obs. 1 




o 


e obs. 2 




□ 


e obs. 3 




A 


^obs. 4 




"avg. 



0.2 
0.0 



x 



o 



A 



-ra- 



1 10 

Vulcanoid Diameter (km) for p R =0.05 



100 



Fig. 6. — Efficiency of detecting synthetic Vulcanoids as a function of object diameter. The 
synthetic Vulcanoids are assumed to have an R-band albedo of 0.05 and a Mercury-like 
spectrum and phase function. For objects with a different assumed albedo, the diameter 

— 1/2 

scales as p R . The observed detection efficiency for objects larger than 8 km is 1.0, while 
the observed detection efficiency for objects smaller than 4 km is 0.0. la error bars are 
shown, based on the Monte Carlo analysis of the data in Table [T] described in the text. 
Given the average detection efficiency of 0.8 in the 5.7-8.0 km size bin, the probability of at 
least one observer detecting a real Vulcanoid of this size is 0.9984. 



